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polynomials corresponding to the Pfaffian. As a particular example we consider the 
chiral unitary ensembles in great detail. 
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1. Introduction 

Random matrix ensembles serve as simple models in a wide range of applications 
[U El El [5l E] which can be found in number theory |7l|8j, disordered systems P, quantum 
chaos [9], empirical data analysis [101 El 112], information theory [13], and quantum 
chromodynamics (QCD) [13]. The complexity of most systems prevents derivations 
of correlation functions whereas analytic results are accessible for the corresponding 
random matrix model. The reason for the applicability of random matrix theory lies in 
the universality of spectral statistics on certain scales like the local mean level spacing 
[T5| [T6| [T7] or on the global scale [TH [191 ISQl [21]. If the Lagrangian of the physical 
system drastically simplifies such that it is effectively described by global symmetries 
there might be a random matrix model fulfilling the same symmetries. 

Already in the 60's and 70's [221 [2S1 [231 [23 [2S], the /c-point correlation functions 
of the Gaussian and circular ensembles for the three symmetries of orthogonal (/3 = 1; 
GOE/COE), unitary {/3 = 2; CUE/CUE) and unitary-symplectic (/3 = 4; GSE/CSE) 
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invariance were derived. They can be expressed as a single determinant for the unitary 
case and a single Pfaffian for /3 G {1,4} where the integrals are pulled inside of these 
structures. Their matrix elements only depend on two eigenvalues which is a drastic 
simplification of the integrand. Since then many other random matrix ensembles were 
studied, e.g. the Ginibre ensembles [271 EH |29l [30] and the the other two rotation groups 
O(A^) and USp (2A^) [31]. The fc-point correlation functions as well as the averages 
over ratios of characteristic polynomials for many of these ensembles are determinants 
and Pfaffians with relatively simple entries only depending on one or two eigenvalues 
[S2| [31] • For a long time it was thought that determinants appear for ensembles with 
P = 2 and Pfaffians for the other two cases. In Refs. [351 [36] the general conditions where 
derived to find these structures. Thus all these particular random matrix ensembles were 
unified in one procedure to derive these structures. 

Very recently a random matrix model for the Wilson Dirac operator was introduced 
[S7] in lattice QCD. It generalizes the chiral GUE which was studied in a Hermitian 
version [371 [381 ISHl HQ] and a non-Hermitian one [H]. The eigenvalue correlations exhibit 
Pfaffians for the Hermitian [ID] as well as for the non-Hermitian case [12] reflecting the 
structure found in Ref. [36]. This structure has to be also valid in the continuum limit 
which is the chiral GUE. Hence the question arises if the Pfaffian determinants obtained 
for the fc-point correlation functions and thus for the averages over ratios of characteristic 
polynomials are much more general than conjectured in the broad literature. 

Also in other intermediate random matrix ensembles Pfaffians were found. For 
example a similar situation arises in the transition from GUE to GOE or GSE [l3l [35] . 
If the ensemble is purely a GUE then then the eigenvalue correlations can be cast into 
determinants whereas the smallest interaction with a GOE or a GSE yields a Pfaffian. It 
would be of theoretical, technical and numerical interest if all ensembles corresponding to 
P = 2 exhibit this phenomenon when coupling it to another random matrix ensemble. 
Such a property simplifies the spectral statistics of intermediate ensembles onto the 
behavior of the entries of the Pfaffian which are averages of one or two characteristic 
polynomials only. 

We derive Pfaffian determinants for averages over ratios of characteristic 
polynomials weighted by a joint probability density function factorizing in weights of 
the single eigenvalues apart from a squared Vandermonde determinant. This squared 
Vandermonde determinant can be cast into one determinant similar to the /3 = 4 case. 
Thus it fulfills the same condition as presented in Ref. [^ which implies a Pfaffian. 
This unifies all ten symmetry classes in the Cartan classification [HI US] and exhibits a 
hidden universal algebraic property in all of these ensembles. 

An introduction of the main idea and of the important functions for the technique 
used here is given in Sec. [H In Sec. [31 we recall some basics known about the 
determinantal structure obtained for averages over ratios of characteristic polynomials 
with respect to chiral unitary random matrix ensembles. In contrast to this structure 
we derive Pfaffians for the same correlation functions in Sec. [H Thereby we discuss 
the Wilson-Dirac random matrix ensemble as a neat application and a good motivation 
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of the derived Pfaffian determinant at the end of this section. The skew-orthogonal 
polynomials corresponding to the Pfaffian determinants are indeed closely related to the 
orthogonal polynomials which are found in the determinantal structures. This relation 
is shown in Sec. O In Sec. El we discuss the generalization of these results for chiral 
unitary ensembles to other random matrix ensembles like GUE and CUE. 

2. Preliminaries 

Structures found in supersymmetry are the key ingredient for the technique used in the 
ensuing sections. These structures allow to derive determinants as well as Pfaffians of 
averaged ratios of characteristic polynomials and, thus, fc-point correlation functions for 
a large class of random matrix ensembles in a direct way. The main idea is to recognize 
that these structures are a pure algebraic property of the random matrix ensemble 
and not an analytic one. By an algebraic rearrangement of the integrand one gets the 
determinants and Pfaffians without explicitly calculating any integrals. This idea was 
first proposed in Refs. [351136] . 

The requirements to obtain determinants was traced back to a factorization of 
the probability density of the random matrix ensemble into densities for the single 
eigenvalues times two Vandermonde determinants (see Ref. [35]), i.e. the measure for 
the single eigenvalues has to be 

N 

d^^{z) = \[g,{z,)\^^{z)\' (2.1) 

i=i 

with the Vandermonde determinant 

A^(.) = W {z^- z,) = (-1)^(^-1)/^ det K-i] . (2.2) 

l<a<6<Af 

The variables z can be complex which correspond to ensembles related to biorthogonal 
polynomials [16]. For Pfaffians this requirement changes to a weight for pairs of 
eigenvalues and a single Vandermonde determinant [36], i.e. 

N 

dfl{z) = Y\_92{z2j-1, Z2j)A2N{z). (2.3) 

i=i 

If one of these two conditions are fulfilled then the technique presented in Refs. [351 136] 
circumvents the integration theorem by Dyson and Metha [251 ESI El HZ] • Moreover the 
approach of Refs. [351 ES] makes an integration theorem unnecessary at the end since it 
is automatically fulfilled for random matrix ensembles traced back to measures of the 
form (12.11) or (12. 3p . This can be readily seen by the combination of the determinantal 
and Pfaffian factorization for averages over ratios of characteristic polynomials [S51I55] . 
the representation of the orthogonal and skew-orthogonal polynomials as averages of the 
corresponding ensemble [HI HHl SHI El SZl ED] and the expressions of the kernels of the 
determinants and Pfaffians in orthogonal and skew-orthogonal polynomials [SI HZ]. In 
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Sections [3] and m we derive the fc-point correlation function without using the integration 
theorem by Dyson and Metha. 

Although, we do not explicitly need supersymmetry, in particular a superspace, 
some functions are quite useful to write the algebraic expressions of the calculations in 
a very compact, constructive and intuitive way. These functions have their origin in the 
theory of supermatrices. For the interested reader, good introductions in supersymmetry 
are given in Ref. [SI] and in the appendix of Ref. [1]. Here we only recall some of these 
useful algebraic functions and notions. 

A diagonal {pjq) x {p/q) supermatrix x consists of two blocks, x = diag(xi,X2). 
The p X p matrix xi and the q x q matrix X2 are indeed diagonal, too. The supertrace 

Str and the superdeterminant Sdet of s is then defined by 

p 1 

Str a; = trxi — tra;2 = ^^Xji — ^^Xj2, (2.4) 



Sdet X 



det xi 



re 



=1 "''ii 



det X2 riLi ' 
The crucial function of the method used here is 

Ap(xi)Ag(x2) 



B 



p/q\ 



a,b 



Xb2) 



(2.5) 



'_l)9('?-i)/2+(g+i)p (jg^ 



for p < q. It is the square root of a Berezinian, 



Xal ~ Xb2 



{^b2 } 



l<a<p 
l<fe<ij 



l<a<g— p 
l<b<5 



Bj/,(x) 



Ber^/J^(x) 



(2.6) 

which is the Jacobian in superspace when diagonalizing a Hermitian (p/q) x {p/q) 
supermatrix. The notation on the right hand side of Eq. (12. 6p refers to the one used in 
Refs. 

Everything we need for the method of Refs. [351 [36] are the functions Sdet and 
B embedded in an ordinary space like WP^'^ or C^"*"^. Hence those readers who are 
not accustomed to supersymmetry may consider these functions as ordinary, rational 
functions. 



3. Review of chiral unitary random matrices 



D 



We consider the ant i- Hermitian random matrix 

W 
-W^ 

which is distributed by the density 

P{D)d[D] = exp[-atTV{WW^)]Y[dReWab d ImWab 

a,b 



(3.1) 
(3.2) 
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with a non-zero normalization constant. In particular it serves as a model for the Dirac 
operator in QCD [13]. The constant a is proportional to n. The matrix is a rax (n + z/) 
rectangular matrix. Each of the n(n + z/) entries of is a complex number which might 
be statistically coupled by the arbitrary density P. The parameter u with < z/ < n is 
the topological charge or also known as index of the Dirac operator such that D has u 
generic zero eigenmodes. The potential V is invariant under the group U (n), i.e. 

ViUWW^U^) = UV{WW^)U\ (3.3) 

and is chosen such that all moments of the ensemble over exist. In the simplest 

case P is Gaussian. Nevertheless the arguments given here are also true for an arbitrary 
potential. We only need the property 

exp[-atTV{A^)] = Y[exp[-aV{X^j)] (3.4) 



P 



\ 






A 


" 


-A 


















for the matrix A = diag (Ai, . . . , A„) with the singular values < Ai < . . . < A„ of W, 
i.e. there are f/ G U (n) and V & U (n + u) with 








A 





D = diag {U, V) 


-A 




















diag(f/^0). 



(3.5) 



In this basis the measure (13. 2p can be written as 

Vol „ Vol 



P{D)d[D] 



Vol ^ Vol 

X dfiTj („)/u "(1) {U)dn\j (n+i/)/u (i/) {y) ■ 



(3.6) 



The abbreviation of the constant 



Vol, 



n 



271^ 



(3.7) 



refers to the volume of the unitary group U(/). Thus, the prefactor in Eq. (13. 6p is 
the volume of the coset [U (n) x U (n + z/)]/[U"(l) x U (z/)]. The measure rf/i© is the 
normalized Haar measure of the coset 0. 

An important quantity to analyze the eigenvalue statistics of this ensemble is the 
average over ratios of characteristic polynomials with respect to D, i.e. 



k2 



n det(L' - iKj2t 



2n+u) 



fci 

Cnx'Cn + i') n det(L) — lHjlt2n+p 



-P{D)d[D] 



(3.1 



with the diagonal, non-degenerate (fci//c2) x (^1/^2) supermatrix k = diag(Ki,K2) = 
diag (kii, . . . , Kfc^i, K12, . . . , Kj!c22) and the 2n + u dimensional unit matrix l2n+u- This 
average is also known as the partition function with ki bosonic and k2 fermionic flavors 
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in QCD [521 ESI [IS]. The variables Kji are complex numbers with a non- vanishing 
imaginary part such that the integral is well defined. The partition function (13.81) 
is simply related to the matrix Green function and, thus, to the /c-point correlation 
function by derivatives with respect to k. 

and 



The joint probability density (13. 6 p is of the class studied in Ref. [35] 
can, therefore, be written as a determinant. This was derived in many articles 
before [SH The crucial idea presented in Ref. ^35] is the combination of the 

ratio of characteristic polynomials (13. 8 p with the two Vandermonde determinants (13.60 
to square roots of Berezinians (12.51) . i.e. 



k2 



n det(A2 - k]^1, 



2i+n\ 



"11 



A2)B 



.A' 



fci 

n det(A2 



K 



Biii//2i('^i)B 



h2/l22\'^lj 



(3.9) 



for any choice of natural numbers lu + /12 = ki and I21 + I22 = ^2- 

In Eq. (13.91) . we split the supermatrix k into the two sets ki = 
diag(Kii,K2i) = diag(fi;ii,...,fi;;iii,Ki2,...,K/2i2) and K2 = diag (ki2, K22) = 
diag (k^j^+i^i, . . . , Kfc^i, K/21+1,2, • • • , /^fc22)- The choice how we split this set is arbitrary 
and, thus, we get equivalent but not trivially related results. This was already rec- 
ognized by the authors of Ref. [55] for products of characteristic polynomials. Let 
di = n + I21 — hi and d2 = n + I22 — /i2- The interesting case is di,d2 > because we 
want to discuss the limit n — )■ 00 and ki, /c2 fixed, at the end of this section. The other 
cases are discussed in Ref. [55] . 

Without loss of generality we assume di < d2- We rearrange the integrand (13.80 
with the help of Eq. (13. 9p which yields 



7(71,1/) 

^ki/k2 



k) oc Sdet 



oc Sdet 



X det 



K 



K 



B 



hi/l2i+n[l-^l, 



A2)B 



h2/l22+n[l^2i 



,A2 



B 



hi/h 



h2/l2 



]^exp[-«r(Af)]Af+MA, 



n exp[-«V(A|)]Af 



BiilA2l('«l 



)B 



h2/l22['^l 



f 1 1 




\ 1 1 


I l^al " H2 , 


l<a<hi 





2(a-l) 
62 



1<6</21 

l<a<di 
l<b</2l 



A 



2{a-l) 
62 



l<a<Zii 
l<6<n 

l<a<di 
l<6<n 



X det 



'^al ^^62 



2(a-l) 
62 



iii+l<a<fci 
i2i+l<6<A:2 

l<a<d2 
«21+1<6<A:2 



~ -^62 

2(a-l) 
62 



/ii + l<a<fci 
l<6<n 

l<a<d2 
l<b<n 



(3.10) 
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Applying the generalized Andreief integration theorem [56l 135] we obtain 
Sdet '"k 



Biil/Z2l(«^l)Bii2/i22('«l) 





(3.11) 



X det 



1 



t^2 ^2 



'2i+l<a<fc2 
«ii+l<fe<fci 



1 



,.2 



K 



al 



^62 



2(a-l) 
b2 



l<a<«ii 
l<f)<i2l 

l<a<di 
l<b</2l 



l<a<hi 



{Fa{>^bl)] 



l<a<di 
iii+l<fe<fcl 



{Ffe(«:ai)} 
[Mab] 



+I<a<fc2 
l<6<(i2 

l<a<Zii 
l<fe<d2 



l<a<di 
l<b<d2 



Notice that Andreief's integration theorem as well as its generalization is only an 
algebraic rearrangement of the integrals without explicitly calculating any integral. The 
functions F and Fa are one dimensional integrals and their explicit expressions are 
not so important as we will see in the discussion below Eq. f l3.15p . For the interested 
reader we refer to Ref. [32] where the explicit integrals are given for general random 
matrix ensembles corresponding to determinants (/3 = 2). The constant di x d2 matrix 
M = [Mab] is given by 



ab 



y2(a+b-2) exp[-«l^(A')]A2'^+MA 



(3.12) 



and thus generates the moments of the measure. 
In the next step we use the identity 

A B 
C D 



det 



AetDdet[A- BD-'C] 



(3.13) 



for arbitrary matrices A, B and C and an invertible matrix D. For the matrix D we 
choose the di x di matrix 

^ = [^^^^]i< b<, (3.14) 

l<a,b<di 

which is only a part of the full rectangular matrix M appearing in Eq. f l3.1ip . The 
determinant of D is proportional to the normalization constant of the ensemble (14. 2 p 
and is therefore invertible. Employing Eq. (I3.13P we find 



</?2(-) 



B«ii/i2l('«l)Bii2A22('«2) 



(3.15) 



X det 



[I^a2, I^b2 



G2 ^\l^al, f^b2] 



'] I 

J l<a<«2i I 



G^2^\K.bl-, f^a2) 



l21+l<b<k2 

/ii+l<a<fci 
l2i+l<b<k2 



G^^'^\Kal, l^bl) 



H'f\nb2] 



di+l<a<d2 
l2l+l<b<k2 



l<a<l2i 
l<b<hi 

/il+l<a<fci 
l<6</ii 



di+l<a<d2 
l<b<hi 



In the last step we identify the functions G^f^\ Gf^\ G^^^\ h["'^ and 112"'^ by considering 
the particular choices (/n, /12, /21, ^22) G {(0, 0, 1, 1), (1, 0, 1, 0), (1, 1, 0, 0), (0, 0, 0, 1), 
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(1, 0, 0, 0)}. In all of these cases the determinant reduces to one of the entries. Then we 
obtain 

di-l 



^0/0 



(_l)fci(fei-l)/2+a2i+l)(A:i+l)+aii+l)(A:2+l) J 



BiilA2l('^l)B«i2A2 



n-1 
j=0 



X det 



^0/2 \'^a27 I^b2] 

"-di-l^O/O 
1 Zfi^^^\Kal, Ki,2) 



1 



(3.16) 



l/l' ('^61) ^02) 



^0/0 ^ "1 



62^ 



^0/1 ^'"^('^bsj 
y(a— 1,1/) 
^0/0 



l<a<«2i 
«2i+l<6<fc2 



iii+l<a<fci 
Z21+l<fe<fc2 



h 



Adi+l,u) ^2/0 
'0/0 



l<a<Z2i 
l<fe<ill 



iii + l<a<A:i 
l<b<ili 



h 



ii+l<a<d2 
l21+l<b<k2 



^(a,u) ^1/0 
^0/0 



cii + l<a<d2 
l<fe<«ii 



for the partition function (13. 8 p which is a particular result of the general one derived in 
Ref. [55] . 

The determinant (13.160 interpolates between one-point and two-point kernels as the 
entries of the determinant. We emphasize again the choice of the numbers < lu < ki 
and < /21 < ki and the splitting of n are arbitrary. The particular choice In = ki and 
hi = yields the ki + ^2 dimensional determinant with one-point kernels considered 
in Refs. [521 [IS]- This choice is suitable for the microscopic limit in chiral random 
matrix theory. For bulk and soft edge correlations [16] the representation in two point 
correlations are the better choice to make contact with other random matrix ensembles 
[T8| fT9| EHl |T7] . This case relates to the choice lu = ki and /21 = k2 for k2 < ki and 
III = and I21 = for k2>ki. 

The fc-point correlation function at the k variables x = diag {xi, . . . , Xk) is given by 

„ n—k 

Rt"\x) oc / A^(diag(x2,A2))exp[-atrF(x2) -atry(A2)]det''^+^x^[Af+MA,• 
JIR^-'= -^^ 

oc A^(a;^)det X exp[— a; tr \/(a;^)]ZQ"2fc'''''^(diag (x, — x)). (3-17) 

Now we employ the formula (I3.16P for (/n, /12, /21, ^22) = (0, 0, k, k) and find the result 



Rt"\x) oc det 



V^exp[-a(V(x2) + V{xl))/2]Z, 



{n-l,u) 
0/2 ' 



Xn 



-Xb 



l<a,b<k 



(3.18) 



Since Z, 



0/2 



{—lyZ^J^ ''^{xa,Xb) this agrees with the general formula for 
(3 = 2 ensembles [5]. Please notice that we derived this formula without using the 
integration theorem by Dyson and Mehta El HZ]. 

The constant h^^'^ in Eq. (I3.16P is the normalization constant of the orthogonal 
polynomial 



i-iy 

^ ^ ^0/0 



(3.19) 
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P^"\x^)pI''\x^)x^''+^ exp[-aV{x^)]dx = h^-^6,i. (3.20) 



The authors of Ref. [57] have shown that these polynomials fulfill a recursion relation 
with respect to the topological charge v by 

p^;^'\x) _ ipt\Q)p%{x)-p%{Q)p^;'\x) 



p(^+i)(0) a:pH(o)pH;(0) -pj.i(0)pf)'(0) 

which is quite useful by taking the limit n — )■ oo. This relation follows when setting 
m = in Eq. (12) of Ref. [57j. One can readily prove identity fl3.2ip by showing the 
orthogonality relation (13.201) for the right hand side with respect to the v + 1 measure, 
i.e. 

Ix 



(3.21) 



x'^ x"^ 





oo 



OC 



E ^^%|^(/>!^^(0)j>a(x^) - pa(0)p;^)(x^))x-^^e-^(^^)rfx 



a=l "-o- 





OC 6,1, (3.22) 

where we used the Christoffel-Darboux formula. The monic normalization of p^j\x) = 
x^ + . . . for all j and u explains the choice of the constants. 

The Cauchy transform of p^^^ is related to the partition function with one bosonic 
fiavor by 

ff;\x')= J ^^^X^'^^^exp[-aViX')]dX (3.23) 



^0/0 

In the result (13.161) we recognize that the choices {In, I12, 121, ^22) = (0, 0, 1, 1), (1, 0, 1, 0), 
(1,1,0,0) yield the same partition functions as the choices (/n, ^12, ^21, ^22) = 
(0, 0, 0, 2), (1, 0, 0, 1), (2, 0, 0, 0), respectively. Therefore the two-fiavor partition 
functions in Eq. (I3.16P can also be expressed in the orthogonal polynomials (I3.19P and 
their Cauchy transforms (I3.23p . i.e. 

^0/2 ^'''\l^a2, K.b2) _ {-Ka2 Kb2Y 



^0/0 
^(di+l,v) 



^(di-l,u) ^ 



a2 '^b2 



pTu<2) piu^L) 

piH<2) P'd^<2) 



(3.24) 



2/0 ' («al,«fel) 1 1 

— det 



(3.25) 
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^1/1 \l^al-il^b2) 1 



h 



det 



Pdi-li'^ll) Pdi-li'^b2) 



il yH2> 
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(3.26) 



^0/0 "-di-l 

These three relations are aheady well known [51 HZ] . They can also be derived with help 
of the Christoffel-Darboux formula. 

The structure f l3.16p is a general property of ensembles with a joint probability 
density including a squared Vandermonde as considered in Sec. 4.2 of Ref. [35] whereas 
the relations f l3.24p - fl3.26p have to be slightly modified for other ensembles. 

In the microscopic limit the authors of Refs. [THl US] have shown that for a generic 
potential V the orthogonal polynomials and their Cauchy transforms become 

Jv{x) 





( X2 






f 


\ [cnY 



n>l 
(X 



oc X Ky(x) 



(3.27) 
(3.28) 



where c is a constant depending on the potential V . The functions Jy and are the 
Bessel function of the first kind and the modified one of the second kind, respectively. 
Hence in the microscopic limit the partition function (13. 80 is 



y{n,v) f ^ \ ">1 
^fcl/fc2 



cn 



1 



oc 



B«ii//2l('«l)Bii2A22('«2) 

{ll'\Ka2,K,,)} 



(3.29) 



X det 



l<a</2l 1- ^ 'J i<a<l2i 

l2i+l<b<k2 l<b<lii 



(«:,!, /tbs)} {li^\l^al, l^bl)} ^ 

(ii+l<a<fei (ii 



{'^b2Ju+a{'^b2)} 



l21+l<b<k2 

0<a<d2-di-l 
l2i+l<b<k2 



{n^-^Ky+ai'^bl)} 



+l<a<fci 
l<b<hi 



0<a<d2-di-l 
l<6<Zil 



where 

lj)^\fia2, fib2) 
li^\nal, K,b2) 



'ia2Ju-lil^a2)Juil^b2) — l^b2Juil^a2)Ju-lil^b2) 



t--2 i~'2 

'^a2 '^b2 

Jv+l{l^a2)Jv-l{l^a2) - Jl{l^a2) 
2 

l^alKy-i{Kal)Jy{Ki,2) - l^b2Ky{hial) Ju-l{.l^b2) 



a ^ h, 
a = b. 



62 



bl 



a^b, 
a = b. 



(3.30) 
(3.31) 
(3.32) 



This is the well known result found in the literature 



4. Derivation of the PfafRan determinant 



In subsection 14.11 we derive a Pfaffian determinant for the same class of chiral random 
matrix ensembles discussed in Sec. [31 A neat application of this Pfaffian is presented in 
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subsection 14.21 This example is the random matrix model for the Wilson-Dirac operator 
in lattice QCD [SZl ESI ESI SHI SI] • 

4-1. Pfaffian determinants in chiral random matrix theory 

We show that the representations in determinants fl3.16p are not the only existing ones 
for chiral unitary ensembles. A non-trivial Pfaffian can be derived for the partition 
function by noticing that the square of the Vandermonde in the measure f l3.6p can be 
rewritten as one Vandermonde determinant of the variables ±Aj, i.e. 

Al(A') = (_i)»(»-i)/2 ^2n(A,-A) 
™^ ^ ^ ^ 2"detA ^ ^ 

The determinant of A will be put into the weight later on, cf. Eqs. (14. 2 p and (14. 3 p below. 
Considering the Wilson random matrix theory [371 ESI ESI HSl E] such a splitting arises 
in a natural way for finite lattice spacing. Then an eigenvalue pair ±«Aj becomes either 
a complex conjugated pair or two independent real eigenvalues corresponding to a pair 
of eigenvectors with positive and negative chirality. Hence, the Pfaffian resulting from 
the single Vandermonde determinant (14. ip is the one which is generalized to non-zero 
lattice spacing and not the determinant (JOj |12] . 

This allows us to define an anti-symmetric two-point measure on 

V{xl) + V{xl] 



g{xi,X2) = — - — exp 



-a- 



6ixi+X2)[Qix^)-Qix2)], (4.2) 



where B is the Heaviside distribution. Then we consider the measure 



n 



m = ^V^^^^2n(A) n^?(A2,-l, X2,)dX2,dX2,-, (4.3) 
1 j=l 

over 2n independent eigenvalues instead of the measure (13. 6p . This measure fulfills the 
general condition for finding a Pfaffian, cf. Ref. [SS] and see also Eq. (12.31) . 

The partition function (13. 8 p can be expressed in terms of this new measure, 

(^^n{k^+k2) r ^" ni/^iZ-Aa) 

= Sdet -^(-m) J n f m- (4.4) 

IR2n n {Kjl - Aa) 

i=i 

In the first step we extend the Vandermonde determinant (14. ip with the characteristic 
polynomials. 



X 



B ki/k2+2n{l^, A) 
Bfci/fc2('«) 



n5'(A2j, \2j-i)dX2jdX 



2i-i 



where Sdet is the superdeterminant. This representation is apart from the 2;2Ar+i-integral 
of the form as in Eq. (3.3) in Ref. [56] • Notice that in this extension we do not have the 
same freedom as in the determinantal case (13.91) since there is only one Vandermonde 
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determinant in the integrand (14.31) . Let d = 2n + k2 — ki > 0. Then we employ the 
representation of the function B as a determinant, see Eq. (12. 5p . 

Sdet 



(4.6) 



X det 



r 1 1 




\ 1 1 






\<a<ki 




l<a<ki 




l<b<k2 




l<b<2n 



<d 
l<b<k2 



<d 
l<b<2n 



The generahzed de Bruijn integration theorem [SHIES] can be apphed now which yields 

Sdet ~'^K, 



7{n,v) 
^kxjki ' 



Zl 'n'. IKI OC 



B 



(4.7) 



X Pf 



{-I2'] 



l<a<fci 
\<b<k2 



l<a<d 
l<b<k2 



-FaiUbl) 



l<a<k2 
l<b<ki 



l<a,b<ki 



l<a<d 
l<b<ki 



Mab 



l<a<k2 
l<b<d 



l<a<fci 
l<b<d 

l<a,b<d 



As Andreief's integration theorem the generalized de Bruijn integration theorem is 
only an algebraic rearrangement of the integrals without calculating any of them. The 
functions F and are two-fold integrals and are again not much of importance, see the 
discussion after Eq. (I4.10p . Explicit expressions of them are given in Ref . ^6] for general 
random matrix ensembles corresponding to Pfaffians comprising the measure (14. 3p . too. 
The d X d anti-symmetric matrix M = [Mab] consists of the moments 



ab 



(Ar^A^^ - Ar^Ar^)(?(Ai, X2)dXidX2. 



Analogously to Eq. (IS.lSp . we employ the identity 



Pf 



A B 

-B^ C 



Pf CPf [A + BC-^B'^' 



(4.8) 



(4.9) 



with an arbitrary matrix S, an arbitrary antisymmetric matrix A and an arbitrary, even 
dimensional, antisymmetric matrix C which has to be invertible. Let ki + k2 be even. 
Then d is also even and the Pfaffian of the matrix M is proportional to the normalization 
constant of the ensemble (13. 2p . Hence the choice C = M is well-defined. This yields 



1 



B 



fci/fcal 



(4.10) 



X Pf 



\l^a2i I^b2] 



-G2 \l^ali f^b2] 



l<a,b<k2 



l<a<ki 
l<b<k2 



G2 \K,bl, I^a2] 



G^3 \l^al, f^bl] 



l<a<k2 
l<b<ki 

l<a,b<ki 
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The functions G^i \ G^2^ and G^^^ can be obtained by considering the cases (^1/^2) = 
(0/2), (1/1), (2/0), respectively. In each of these cases the Pfaffian (14.101) reduces to a 
single term. This leads to a particular case of the general result derived in Ref. [3S]- We 
find our main result of this article 



d-l 



y{n,u) 
^0/0 



n-1 
j=0 



(4.11) 



X Pf 



K'b2 — Ha2 



"■d/2-1^0/0 



^0/2 \l^a2^l^b2j 



7{d/2,u) 
'l/l ' 



^al5 >^b2) 



y{d/2,v) 

"0/0 



[Kal - Kb2 



Ad/2,u) 
'1/1 ' 



l^bl, I^a2) 



y{d/2,v) 

^0/0 



[K.a2 - /«W 



^1/2 ('^al 



t^bli 



y(d/2+l,v) 
^/0 



y{d/2+l,v) 
'2/0 ' 



for even ki + The indices a and 6 run from 1 to k2 in the first columns and the first 
rows and from 1 to ki in the last ones. The result for odd ^2 + fci can be readily obtained 
by introducing an additional fermionic flavor and sending it to infinity. This shifts the 
parameter ci to c? + 1 and adds a row and a column to the matrix in the Pfaffian (14. lip 



with the partition functions Z. 



{{d-l)/2,u)^ 
0/1 



Kb2) and Zj^g"^^^^^''''' 



(ubi) which are apart from 



a factor k'^ an orthogonal polynomial and its Cauchy-transform, cf. Eqs. fl3.19p and 
(I3.23p . Notice that the matrix in the Pfaffian (14. lip is indeed antisymmetric because 
'^ti2~^'^^ and Z^^j^'^^''^^ are symmetric under a permutation of the entries. 
Indeed, Eq. (14.111) cannot be traced back to the identity 



Pf 



X 

-X^ 



'_^Wp-l)/2 ^g^^ 



(4.12) 



with an arbitrary p x p matrix X. We refer to the relation (I4.12p as a trivial Pfaffian 
extension of a determinant. The Pfaffian (14. lip seems to be the result of recursion 
relations of the orthogonal polynomials (I3.19p . It is difficult to see how these recursions 
have to be performed to map the Pfaffian (14.110 to the determinant (I3.16p . However the 
construction of this structure seems to be the same for a broad class of ensembles. This is 
confirmed by the fact that the result (14. lip can be extended to all factorizing ensembles 
with a squared Vandermonde determinant in the joint probability density (12. ip . This 
will be shown in Sec. El 

Again one can consider the A;-point correlation function (I3.17P and what it looks like 
with the Pfaffian determinant. Using the result (14. lip we find for the /c-point correlation 
function 



R 



X 



oc exp[— a tr V"(x)] 



X Pf 



/ \ y(n—l,u) / \ 



1 Xn 



Xb)Z. 



Jn^^l^. 
0/2 ' 



-Xb) 



Xb)Z^ 



{n-l,u}, 
0/2 ^ 



-Xb) 



l<a,b<k 



(X exp[— a tr V{x)] 
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(4.13) 



l<a,b<fc 



Again we have not employed the integration theorem by Dyson and Mehta [251 [26| l5| H7] . 
To see that Eq. fl4.13p indeed agrees with the determinant (13 .181) one can consider the 
square of the Pfaffian, 



Pf ' 



det 



2Met 



2^^ det 



22'=det2xdet^ 



(X(i X\)')Zi^j2 (-^aj -^b) 


(X(j -|- Xf^ Z Qi2 (3^05 •^b) 


^"^a ~^ ^b)Zf^l2 {Xa, Xh) 


( \ yl""^'*^) ( \ 
\Xa •^b)^Qj2 y^ai^b) 

(X(j -|- X\^Z Q 12 (^^a; -^b) 







l<a,6<fc 



l<a,fe<fc 



0/2 



[Xa + 



0/2 ' 



0/2 



\Xa X}yjZjQi2 \Xa)Xl)j 



J l<a,6<fc 



7(n-l,i/) 
■^''^0/2 











(^a ~l" '^b)-^o/2 








l<a,6<fc 



^0/2 (a-a) -^fc) 



J l<a,fe<fe 

The square root of Eq. ( Kl^ yields Eq. f lXTSj) . 

In the large n limit, we employ Eqs. (l3.24tlX^S]l and (13.30113.3^ and obtain 

(n,u) / ^ ">1 1 



(4.14) 



era 



B 



ki/k-z 



(4.15) 



X Pf 

for even ki + ^2 and 



[K,a2 — f^b2 



I^a2, I^b2) 



[f^al + fib2)li^\fial,fib2] 



[l^bl + f^a2)I^\l^bl, I^a2) 



al 



l^bl 



ki/k2 



cn 



ri>l 

oc 



B 



ki/k-z 



(4.16) 



X Pf 



-Jv{l^a2) 



-KJKal] 



Ju{l^b2) 



I^a2, I^b2j 



{Kbl + Ka2)ll^\ 



l^bl, I^a2] 



-{Kal + I^b2)li^\l^al, «b2) 

for odd ki + k2- These Pfaffians carry over to the Wilson Dirac random matrix 
model lini II2] . For small numbers of bosonic and fermionic flavors these results were 
checked by the recursion relations of the Bessel functions [59J. 

Please notice the difference in the prefactor of Eqs. (I3.29p . (I4.15P and (I4.16p . 
The entries of the Berezinian are the squares of the variables k for the determinantal 
structure (13.291) whereas it is only k, for the Pfaffian. This yields a technical advantage 
when calculating eigenvalue correlations of the random matrix models for the Wilson 
Dirac operator. 
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4-.2. An application: Wilson- Dirac random matrix 



The Wilson-Dirac operator is a modified Dirac operator on a lattice. In the infrared 
limit this operator can be modeled by the Wilson-Dirac random matrix [371 EH ESI HI] 
which is a {2n + u) x (2n + u) Hermitian matrix 

aA W 
-W^ aB 



distributed by the Gaussian 



(4.17) 



P{D 



wj 



exp 



n 



-{tr A^ + tr B^) - ntiWW^ 



(4.18) 



The variable a plays the role of the lattice spacing. The chiral symmetry is explicitly 
broken by the Hermitian matrices A and B, i.e. 

75 ^wL=o75 7^ - ^wL=o ^i*'^ 75 = diag(l„, -!„+!,), (4.19) 
which have the dimensions n x n and {n + u) x [n + u), respectively. Hence, A and B 
model the Wilson-term. 

We consider the partition function with Nf fermionic flavors, 

Nf 



m, a) = / n det(£'w + mjl2n+u)P{D^N)d[D.f^]. (4.20) 

The external variables m = diag (mi, . . . , m^^) play the role of the quark masses. Indeed 
one can also consider bosonic flavors. However, we restrict ourself to fermionic flavors 
to keep the example as simple as possible 

In the microscopic limit (n — )• oo), m = 2nm, a = y/na/2 and u are kept fixed. 
This yields the integral 



exp 



- tr m{U + U^) - tr [U^ + {U^f) 



defUdfiiU). (4.21) 



V{Nf) 



For a derivation of this result we refer to Refs. |37[ EH]- Exactly the integral fl4.21|] 
makes contact with lattice QCD [60l EH [62], 163] . 

At zero lattice spacing (a = 0) this partition function can be identified with the 
one considered in Sec. [31 



m,a = 0) Qc 



O/Nf ' 



tm). 



(4.22) 



Considering again the microscopic limit fl4.2ip . we trace the integral back to the a = 
result by introducing a N{ x Nf Hermitian random matrix a similar to the calculation 
in Ref. [23151], 

1 



Z]v^ '{m,a) oc 



exp 



4a2 



trier — im) 



2{aNi 



X 



oc 



exp [-ttr a{U + U^)] det'' Udfi{U)d[a] 



\J{Nf) 



exp 



4a2 



tr(cr + tm 



)2 - 2(aiVf)2 



Zi7j;^,{a)d[a]. 



(4.23) 
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Notice that a is an ordinary matrix and not a supermatrix because we consider fermionic 
flavors, only. The constant exp[— 2(aA^f)^] can be shifted into the normalization constant 
and can, thus, be omitted in the ensuing calculations. 

A diagonalization of a = UsW with G U (A^f) yields a Harish-Chandra-Itzykson- 
Zuber-integral [651 166] in the Gaussian term. The partition function zj^^'^^ is invariant 
under U(iVf). We find 



Z]v^ '{m,a) oc 



det [exp [-{sj - tmif/4a% 



Atvj (m) 



(4.24) 



Employing the result as a determinant of the microscopic limit of a = partition 
function, cf. Eq. (13. 29^ . we end up with a complicated expression, 

det [exp \-{sj - innf / Aa^]]^^..^^^ 



Z]vj '(m, a) oc 



n>l 



det [srV;.-i+j(s)]i 



Ajv,(s) 



■d[s\ 



(4.25) 



There is no obvious way to further simplify the integral (14.251) due to the factor 
AjVj(s)/AjVj(s^). This was not much of a problem for the authors of Refs. El] 
because they only considered a small numbers of flavors. However the problem is highly 
non-trivial for an arbitrary number of flavors. 

This problem can be solved by using the Pfaffian expressions (I4.15P and (14.161) 
instead of the determinant. Let A^f be even to keep the expressions as simple as possible. 
Then we have for the microscopic limit (I4.2ip 



Z]v^ '[m,a) oc 

X Pf 



n>l 



det [exp [— (sj — ifhiY /Aa^]]^ 



<j,i<N{ 



A7Vj.(m) 

Sj ~h Si 



(4.26) 



d[s\ 



J l<j,i<Nt 



After expanding the determinant no term hinders us to pull the integrals into the 
Pfaffian. We obtain the compact result 

zl;-^(m,a)oc^-i^Pf 



with 

Z2 [mi, m2, a) oc 



mi — 777,2 



exp 



[si - ifhif + (s2 - 7777,2)^"' 
4^2 



(4.27) 



X ^ — - — ^ — ^ — ^ — -dsids2. 

Si + S2 

This is a drastic simplification of the problem compared to Eq. (I4.25p . 



5. Skew-orthogonal polynomials 



What are the skew-orthogonal polynomials which correspond to the Pfaffian (14. lip ? 
In order to solve this problem we consider the two-point measure (14. 2p . The skew 
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orthogonal polynomials qj are defined by 

q2i-i{xi) q2i-i{x2) 



det 



g{xi,X2)dxidx2 



det 



q2'j{xi) q2j{x2) 
q2i{xi) q2i{x2) 



g{xi,X2)dxidx2 = 0, 



and 



det 



(l2j+l{Xl) q2j+i{x2) 
q2i{xi) q2i{x2) 



g{xi,X2)dxidx2 = hf'^dij. 



(5.1) 



(5.2) 



Moreover one has to assume that qi is a polynomial of order /. 

The integral over the measure (14. 2 p for two arbitrary and conveniently integrable 
functions /i and /2 can be simplified to 

/l(xi) fi{x2) 



det 



det 



/2(a;i) /2(a;2) 
Hx) M-x) 

f2{x) f2{-x) 



g{xi, X2)dxidx2 



X 



2u 



exp [— nF(a;^)] dx. 



(5.3) 



Due to this identity the skew-orthogonal polynomials q['^^ are related by the orthogonal 
polynomials pi in the following way 

qti'{x)=pt\x') (5.4) 
for the even polynomials and 

q^i\i{x) = xp\'^\x'^) + const. p\^\x'^) (5.5) 

for the odd polynomials. Notice that these skew-orthogonal polynomials for V{x) = x 
(the Laguerre ensemble) are similar to but not completely the same as the one for /3 = 1 
and /3 = 4 shown in Ref. [5l HT] for the Laguerre ensemble. The reason is the two point 
weight which is 

Xi — X2 



fl'chGOE(a;i,a;2) = (xia;2)''exp [-a{xl + xl)] 



(5.6) 
(5.7) 



\Xi — X2\ 

fi'chGSE(a;i,a;2) = (xiX2)^''"*'^/^ exp [-a{xl + xl)] 6\xi -X2), 

in comparison see Eq. fl4.2p for /3 = 2. The labels "chGOE" and "chGSE" refer to 
the chiral Gaussian orthogonal ensemble {(3 = 1) and to the chiral Gaussian symplectic 
ensemble (/3 = 4), respectively. The sign function (xi — X2)/|xi — X2I generate the 
modulus of the Vandermonde determinant for (3 = 1. The distribution 6' is the first 
derivative of the Dirac delta function and cancels with these terms of the Vandermonde 
determinant which are zero at the support of the Dirac delta functions. This generates 
Cramers degeneracy in the quaternion case (/3 = 4). 

The solution of Eqs. (15. ip and (15. 2p is not unique which is reflected by the arbitrary 
constant in the odd polynomials (15. 5p . One can readily confirm that this choice of 
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the polynomials solves the conditions (15. ip and (15. 2p by recognizing the symmetry 
qj{—x) = {~iyqj{x) and the orthogonality relation (I3.20p for pj. The normalization 
constant is 

h^^ = hf\ (5.8) 

This relation between orthogonal and skew-orthogonal polynomials seems so trivial 
because of the particular and simple structure of the two-point weight (14. 2p . 



6. A few more ensembles with Dyson index /3 = 2 and PfafRans 

The algebraic rearrangement for chiral unitary ensembles described in Sec. H] can 
be extended to other random matrix ensembles which have a squared Vandermonde 
determinant in the joint probability density function. By the same trick as in Eq. (14. ip 
we write 

Al{z) = (-l)"^"-^^/'^^fl;^=^' (6-1) 

where the variables z = diag^zi, . . . , zn) might be complex. The square root is the 
positive one but this is without loss of generality since the right hand side of Eq. (16. ip 
comprises both roots. Again the determinant of z will be put to the measure dfi for a 
single eigenvalue. 

We consider an average over ratios of characteristic polynomials for random 
ensembles like GUE and CUE, i.e. 

fc2 

N IK^^i-M 

^SU'^) = J ^Ni^) n '-^ — ^/^(^^)' (6-2) 

where dfi is a measure on C and k is chosen such that the integrals exist. Notice that 
there is no modulus of the Vandermonde determinant which is a necessary property of 
the following discussion. A modulus of the Vandermonde is an obstacle to map Eq. (16.21) 
to the general joint probability density corresponding to the Pfaffian, see Ref. |36], which 
we have not managed yet. A modulus corresponds to the biorthogonal polynomials [IS] 
whereas the choice without the modulus corresponds to the orthogonal polynomials, 
only. Apart from the modulus of the Vandermonde it is exactly the correlation function 
discussed in Sec. 4.2 of Ref. |35j. 

With the help of the derivation in Sec. H] the integral (16. 2p can be written as 

- b;:;^ («-^' 
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X Pf 



^0/0 



1 Zy\{Kal,Kh2) 
^0/0 ^ V / 



1 /ta2) 



^0/0 ^ ^ 



for fci + k2 even and 



X Pf 







(6.4) 




7{d) 
'0/0 



1/0 ' 



d 



^0/0 



0/2 \f^a2,f^b2) 



'0/0 



^/O 



35^^J/o('«6l) 
"0/0 



1 -2^^/1(^61,^02) 
(y^ - 

^0/0 ^ ^ 



^^al - yft^bl 7^{d+l)( X 
(J) ^2/0 I'^alj'^Wy 



(rf+l)Z, 



0/0 



for /ci + ^2 odd. The variable c? is 

N+{k2-ki)/2, A;2 + /£ie2N, 
AT + (A;2 - A;i + l)/2, A;2 + /ci + 1 G 2M. 



(6.5) 



The indices a and 6 run from 1 to ki for Ki and from 1 to /c2 for /€2. Apart from the 
square roots of the variables /t these structures are exactly the same as those of random 
matrix ensembles with Dyson index /3 G {1,4}. Hence, it seems to be that the Pfaffian 
determinants (16. 3p and (16. 4p for the average over characteristic polynomials are more 
general than the determinant derived in Ref. |35j . 

Random matrix ensembles whose generating functions can be cast into the 
form (16. 2 p have this non-trivial expression as a Pfaffian. The Hermitian Gaussian 
unitary ensemble as well as its generalization with other potentials fulfill a priori this 
requirement since the joint probability density has a squared Vandermonde determinant 
without the modulus. More generally our derivation applies to each ensemble with 
a real spectrum, a squared Vandermonde determinant and a factorizing probability 
distribution, cf. Eq. (13. 4p . Also the CUE (unitary group) can be cast into the form (16. 2p . 
More ensembles can be found in the tables 1 and 2 of Ref. [35]. The Ginibre ensemble 
as well as its chiral counterpart are not in this class. Their joint probability density 
incorporates a modulus of the Vandermonde determinant and is, thus, in the class for 
the bi-orthogonal polynomials. Therefore it is possible that their eigenvalue correlation 
functions cannot be expressed in Pfaffians like Eq. (16. 3p and (16.41) . 

The skew-orthogonal polynomials corresponding to the Pfaffians (16. 3p and (16. 4p 
have the same relation to the orthogonal polynomials as chiral unitary ensembles, see 
Eqs. (I5.4tl5.8p . By construction this relation is so simple. 
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7. Remarks and conclusions 

We derived a non-trivial Pfaffian determinant for the average over ratios of characteristic 
polynomials of a large class of random matrix ensembles with Dyson index /3 = 2. This 
structure is similar to the one for /3 G {1,4}, cf. Ref. |36]. Hence, it is universal 
and unifies most of the symmetry classes known in the literature, particularly the 
Cartan classification [HI US]- It is unclear how far beyond this classification [ST] this 
structure is applicable. It is only known that there are some of them which share the 
identity (14. lip . For example the real and quaternion Ginibre ensembles as well as there 
chiral counterpart fulfill an identity similar to Eq. (14. lip . 

For many random matrix ensembles like the GUE it seems an academical question 
if one can derive a Pfaffian or not since there are no applications, yet. However, for the 
chiral GUE it is important to know this due to the new results obtained for the Wilson 
Dirac random matrix ensemble discussed in Refs. [551 ESI SOI El 1121 Ell- Pfaffians 
were found there for finite lattice spacing. On the level of joint probability density the 
authors of Ref. [10] checked that the ensemble is the chiral GUE as well as the GUE at 
certain values of the lattice spacing. However for the eigenvalue correlation functions 
the continuum limit has not yielded the known determinant fl3.29p . With this work we 
clarified this puzzle. 

For intermediate ensembles in general our result might be helpful to understand the 
structure appearing by switching the interaction between the two ensembles on. It is 
numerically advantageous to think about spectral correlations of intermediate ensembles 
as kernels of Pfaffians since the integrand drastically simplifies. In combination with 
the supersymmetry method [HSl ESI EOl E] also the number of integrals reduces a lot. 

Moreover, we derived the relation between the orthogonal polynomials and the 
skew-orthogonal polynomials corresponding to the determinants and the Pfaffians, 
respectively. This relation, see Eqs. (15. 4115. 8p . is not only quite simple but also universal 
since it applies to all random matrix ensembles discussed in this work. The relation 
between the orthogonal and skew-orthogonal polynomials for P = 2 slightly differs to 
those found in Ref. [1^ for the cases /3 = 1 and /3 = 4. The difference in the two-point 
weight is the reason for this. Based on the representations (13.160 . (16. 3p and (16. 4p shared 
by all random matrix ensembles with /3 = 2 as well as checks of these representations 
[5S], we conjecture that the recursion relation of the orthogonal polynomials connects 
the determinant and the Pfaffian and this has to be done in a general way. 

The Pfaffian found for the average over characteristic polynomials carries over to 
the /c-point correlation functions. This structure is valid in the large matrix limit, too. 
It should not depend on which scaling limit is chosen since the Pfaffian is independent 
of the matrix size. Hence, the correlation functions appearing as kernels of this Pfaffian 
have non-trivial recursion relations mapping the determinant to the Pfaffian. 
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